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ABSTRACT 

The effects of spatial correlations of density fluctuations on merger histories of dark 
matter haloes (so-called ^merger trees') are analysed. We compare the mass functions 
of dark haloes derived by a new method for calculating merger trees, that proposed 
by Rodrigues & Thomas (RT), with those given by other methods such as the Block 
model, the Press-Schechter formula and our own formula in which the mass functions 
are analytically expressed in a way that takes into consideration the spatial correla- 
tions. It is found that the mass functions given by the new method are well fit by those 
given by our formula. We believe that new method (RT) naturally and correctly takes 
into account the spatial correlations of the density fluctuations due to a calculated, 
grid-based realisation of the density fluctuations and so is very useful for estimating 
the merger tree accurately in a way that takes into consideration spatial correlations. 

Moreover, by applying our formula, we present an analytic expression which re- 
produces the mass function derived by the Block model. We therefore show clearly 
why and how the mass functions given by the new method and the Block model are 
different from each other. Furthermore, we note that the construction of merger trees 
is sensitive to the criterion of collapse and merging of overlapped haloes in cases in 
which two or more haloes happen to overlap. In fact, it is shown that the mass function 
is very much affected when the criterion of overlapping is changed. 

Key words: galaxies: formation - galaxies: mass function - large-scale structure of 
the universe 



1 INTRODUCTION 

Recent observations by the HST and the Keck telescope have 
rapidly increased the available numbers of observations of 
faint, high redshift galaxies, providing us with significant 
new information about the birth and evolution of galaxies. 
In order to understand the significance of these observations, 
it is very important to theoretically understand the physical 
processes underlying galaxy formation. 

In our standard understanding, it is considered that 
dark matter dominates in our Universe and that galaxies and 
clusters of galaxies have formed by the gravitational growth 
of initial small density fluctuations. The fluctuations of dark 
matter have collapsed and virialised by self-gravitational in- 
stability into objects which are called 'dark matter haloes' 
or 'dark haloes'. The larger haloes are generally considered 
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to have formed hierarchically by clustering of smaller haloes 
(it is so called 'hierarchical clustering'). 

The baryonic gas also has collapsed and virialised fol- 
lowing the collapse of the dark matter. Furthermore, in 
the process of galaxy formation, baryonic gas must dissi- 
pate the internal energy by radiative cooling and shrink 
because galaxies are much denser systems than the viri- 
alised dark haloes and stars must be formed inside the cold 
and dense gaseous systems. In analysing the formation and 
evolution of the galaxies, we must treat various physical 
processes over a large dynamic range from ~ 1 — 10^ Mq 
(star formation, heating processes of gases from supernovae, 
dynamical and chemical evolution of gases, etc.) through 
~ ]^q6-i2j^^ (mergers of galaxies, tidal interactions, etc.) 
to ~ Mq (clusters of galaxies, large scale structure 

of the Universe, etc.). So it is difficult to attack the problem 
of formation and evolution of galaxies in a way that connects 
all of the above complicated physical processes. 

One method of analysing galaxy formation is numeri- 
cal simulations (e.g., Navarro & White 1993, Katz & Gunn 
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1991, Katz 1992, Steinmetz & Muller 1994, Evrard, Sum- 
mers & Davis 1994) which directly pursue the gravitational 
growth of dark matter and thermal processes of gaseous sys- 
tems. The advantage of simulation is that we can trace the 
complicated processes of the systems quantitatively. How- 
ever, it is impossible to deal with the very wide range of 
the mass scales simultaneously in the current limited ability 
of computers. Furthermore, since the CPU time is heavily 
consumed in the simulations, it is difficult to analyse galaxy 
formation statistically by pursuing many samples of the sys- 
tem with different initial conditions and parameters in the 
physical processes. So we can say that the problem of galaxy 
formation has never been resolved completely by numerical 
simulations without any uncertainty. 

On the other hand, there is another approach to solv- 
ing galaxy formation, that is, the semi-analytical approach, 
which has been developing. Some works in this approach are 
briefly reviewed below. The pioneer works in this approach 
are, for example, Rees & Ostriker (1977) and White & Rees 
(1978). Rees & Ostriker remarked upon the importance of 
the dissipational process of baryonic gases through radiative 
cooling. They asserted that a dense object like a galaxy must 
have a cooling time scale TcooI shorter than the dynamical 
time scale Tdyn of the system at the virialised stage because 
the object must dissipate the energy effectively in order to 
become a more dense system and malce stars in the cold and 
dense gases. Therefore, one of the criteria for a virialised ob- 
ject to become a galaxy in the future is TcooI < Tdyn at the 
virialised stage. Then Rees & Ostriker showed in the cooling 
diagram that the range TcooI < Tdyn corresponds to the char- 
acteristic mass range of galaxies, i.e., IO^A^q ^ 10"*^^ Mq. 

White & Rees (1978) showed the luminosity function 
of galaxies by using the Press-Schechter formula (Press & 
Schechter 1974;hereafter PS, see below) in which the mass 
function of dark haloes is analytically estimated. Further- 
more, White & Frenk (1991) extended this approach to con- 
sider the merging process of dark haloes approximately by 
using an extended PS formula (Bower 1991). In this formula 
they take into account the conditional probability of finding 
a region with mass Mi at redshift zi which is also included 
in a region with mass M2 at Z2. They calculate the present 
stellar and gas abundance by tracing the merger process of 
the dark haloes from an initial time till the present time in 
each mass scale of the halo on the assumption that there is 
a relation between the halo mass, the gas cooling rate, the 
accretion rate of gas onto the halo and the star formation 
rate. The extension of the PS formula is also discussed by 
Lacey & Cole (1993). By using the peak formula (Peacock & 
Heavens 1985, Bardeen et al. 1986; hereafter BBKS), Lacey 
& Silk (1991, 1993) also investigate the time scale of the 
collapse of galaxy groups after the galaxies have collapsed 
on the assumption that the mass ratio of groups to galaxies 
is constant. 

However, the status of the gaseous systems and for- 
mations of stars in each dark halo at an epoch depend on 
the thermal history of the gaseous systems, merging of the 
galaxies and the merging history of the halo. So it is very 
important to know the merging history of each dark halo 
at an epoch in order to evaluate the status of the gaseous 
systems and the formation rate of stars in the halo at that 
epoch. KaufFman & White (1993; hereafter KW) constructed 
a merger tree that expresses the merging history of the dark 



haloes by considering progenitors of each halo at every time 
by using a Monte Carlo method of the extended PS formula. 
Kauffman, White & Guiderdoni (1993) then calculated the 
statistical properties of galaxies in their model. Another ap- 
proach to construction of a merger tree is the Block model, 
which was developed by Cole & Kaiser (1988). The Block 
model takes the Monte Carlo procedure described as follows: 
First of all, the density contrast is assigned to the largest 
block with mass Mo(~ 10^*^ M0) which has the variance of 
the density contrast, a'^{Mo). Next, the largest block is di- 
vided into two blocks with the same mass Mi{= Mo/2). The 
additional positive density fluctuation generated by a ran- 
dom Gaussian distribution with variance a^{Mi) — a'^{Mo) 
is assigned to one of the two divided blocks with mass 
Mi(= Mo/2) and a negative fluctuation with the same ab- 
solute value of the amplitude as the positive one is assigned 
to the other divided block. The same procedure is repeated 
for Ml, M2, • • • down to the smallest mass scale under con- 
sideration, thereby constructing the merger tree. Cole et al. 
(1994) also calculated the statistical properties of formation 
and evolution of galaxies in their Block model. 

Incidentally, it has been pointed out that the PS for- 
mula has the following crucial weak points. The PS formula 
is derived as follows: In an Einstein-de Sitter universe, spher- 
ical overdense regions collapse and virialise when their linear 
density contrast reaches Sc = 1.69 (see, e.g., Peebles 1993). 
Then, assuming that the density fluctuations obey a random 
Gaussian distribution and that the collapse of the haloes is 
spherically symmetric, we can get the volume fraction of the 
regions of collapsed objects whose masses are larger than the 
mass M, /(> Sc, M). So, the region of the dark haloes with 
mass scale M is equal to df/dM. Thus they proposed the 
formula for counting the number density of objects of mass 
scale M. However, only overdense regions were considered 
in the analysis. Even if the density contrast smoothed on 
the mass scale A4 is less than 5c, there is the case such that 
the density contrast smoothed on the larger mass scale than 
M is greater than Sc- We must consider this case to count 
exactly the number density of the dark haloes. This prob- 
lem is called the 'cloud-in-cloud' problem. Press & Schechter 
(1974) simply multiplied the number density by a 'fudge fac- 
tor' of 2 without a good reason. 

Peacock & Heavens (1990) and Bond et al. (1991) at- 
tempted to solve the cloud-in-cloud problem by using a peak 
formula. They considered the upcrossing probability that 
a density contrast which is lower than 5c with a smooth- 
ing scale M exceeds 5c for the first time when increasing 
the smoothing scale. In the case that the density field is 
smoothed with the sharp fc-space filter, they found that the 
factor of 2 introduced by PS is correct. The cloud-in-cloud 
problem is also solved for Poisson fluctuations by Epstein 
(1983). 

In the cloud-in-cloud problem, however, we must con- 
sider the spatial correlations of the density fluctuations since 
objects have non-zero size in reality. Yano, Nagashima & 
Gouda (1996; hereafter YNG) analysed the cloud-in-cloud 
problem talcing explicitly into account the spatial correla- 
tion of the density fields. They explicitly introduced the 
two-point correlation function in the mass function by using 
Jedamzik's formula (Jedamzik 1995) in which the number 
density of the collapsed objects is given in the form of the 
integral equation which is different from the PS formula. 
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It is also proved by YNG that the results derived from the 
Jedamzik formula are consistent with those derived from the 
PS formula on the assumption that the collapse is spheri- 
cally symmetric when the spatial correlations are not taken 
into consideration. However, YNG showed that the spatial 
correlations greatly affect the mass function. Therefore, we 
believe that the merger tree is also affected by the spatial 
correlations whose effects have never explicitly been taken 
into account in the KW method and the Block model. So we 
believe it is very important to analyse the effects of spatial 
correlations of the density fields on the merger trees of dark 
haloes. 

Recently, a new approach to construction of merger 
trees has been proposed by Rodrigues & Thomas (1995) (we 
call this the Merging Cell model for convenience throughout 
this paper). In their model, the random Gaussian density 
fluctuation field is realised on spatial grids as in construct- 
ing initial conditions of N-body simulations, so it is expected 
that this model naturally includes information about the 
spatial correlation. Then, by finding the region of the col- 
lapsed cells or blocks with each mass scale whose density 
contrast 5 is Sc, we can construct the merger tree (see §2.1). 
Wc expect that this model is more realistic and useful for 
galaxy formation although spherical collapse and random 
Gaussian density fiuctuations are assumed. 

In this paper, we show the mass functions of dark haloes 
by calculating the Merging Cell model. They are compared 
with the mass functions given by the PS formula and the 
Block model. We then explicitly show that the main origin 
of differences between the mass functions given by the dif- 
ferent methods result from the effect of spatial correlations 
by comparing with the mass functions derived by YNG's 
formula (hereafter, the YNG formula). We find that the 
mass functions given by the YNG formula in taking explic- 
itly into account the spatial correlations are consistent with 
those given by the Merging Cell model and so the Merging 
Cell model correctly includes the effects of the spatial cor- 
relations. We believe that this effect is very important for 
calculating the merger tree of the dark haloes. Furthermore, 
by applying the Jedamzik formula, we present an analyti- 
cal expression of the mass function derived from the Block 
model and show quantitatively why and how the mass func- 
tion given by the Block model is different from those derived 
from the PS formula and the Merging Cell model. 

In the Merging Cell model, some haloes happen to over- 
lap since the non-zero size of the haloes is considered explic- 
itly. So, we must consider the criterion of collapse and merg- 
ing for the overlapped haloes. We also show how the criterion 
affects the mass function in the Merging Cell model. 

In §2, the Merging Cell model and the Block model are 
reviewed briefly. In §3, we give the analytical formulae for 
estimating the mass functions by using the Press-Schechter 
formula, the Jedamzik formula and the YNG formula. In §4, 
it is shown that the Merging Cell model is consistent with 
the mass function which is derived by the analytical formula 
in which the spatial correlations are taken into account. The 
mass function given by using the Jedamzik formalism which 
reproduces those given by the Block model is also presented. 
Furthermore, we also show how the mass function is changed 
if the overlapping criterion is changed. We devote §5 to con- 
clusions and discussions. 



2 MODELS OF MERGER TREES 

In this section, we briefly review the Merging Cell model 
(hereafter MCM) and the Block model. 

2.1 Merging Cell model 

We briefly review the MCM according to the procedure and 
the notations shown in Rodrigues & Thomas (1995). 

First, the random Gaussian density field is realised in a 
periodic cubical box of side L. In the random Gaussian dis- 
tribution, Fourier mode of density contrast 5{= {p—p)/p, p is 
density and p is the mean density of the universe) obeys the 
following probability for its amplitude and phase(BBKS), 

P(|^kUk)d|5.|#. = ^ exp |-^| d|5.|^gi, (1) 

where is the random phase of 5k, <5k = |(5k| exp (ii^k) 
and Pik) is the power spectrum (j5k|^), where the angle 
brackets mean the ensemble average of the universe. Then, 
the density contrast at each grid ('cell') is given by Fourier 
transform, 

'5(x) = ^^''5ke*''d^fc, (2) 

where fee is the cut-off wavenumber. 

Next, averaging the density fluctuations within cubical 
blocks of side 2, 4, . . . , L, the fluctuations of the various 
smoothing levels axe constructed. At each smoothing level, 
displacing the smoothing grids by half a blocklength in each 
direction of each axis, eight sets of overlapping grids are 
constructed in order to reflect the position of density peaks 
approximately (see Fig.l). 

Then, the density fluctuations within blocks and cells 
are combined into a single list and ordered in decreasing 
density. The fluctuations are investigated from the top of the 
list. It is decided by the following rules whether each block 
or cell can collapse. Note terminology that halo is a block or 
cell which has already collapsed, and an investigating region 
is a block or cell whose linear density contrast is just equal 
to 5c at the reference time. We investigate whether or not 
an investigating region can collapse at that time according 
to the following rules. 

(a) If an investigating region includes no haloes, the in- 
vestigating region( feZocfe or celt) can collapse and can be iden- 
tified as a new halo. 

(b) When an investigating region includes a part of a halo, 
if the overlapping region has a larger mass than the mini- 
mum of the masses of the halo and the investigating region, 
then the investigating region can collapse (sec Fig. 2). This 
is the criterion of collapse of the investigating region and 
merging for the overlapped haloes. We call this criterion the 
overlapping criterion in this paper. 

(c) In the case of (b), if the investigating region has two 
haloes whose overlapped region within the investigating re- 
gion is more than a half of the region in each halo, the re- 
gion cannot collapse. We set this condition in order to avoid 
long filamentary objects. This is the criterion for linking of 
haloes. 

These criteria are those chosen by Rodrigues & Thomas. 
We will change the overlapping criterion (b) later and see 



© 0000 RAS, MNRAS 000, 000-000 



4 M. Nagashima & N. Gouda 



i+1 level 




Figure 1. Scheme of averaging density fluctuations demonstrated 

in 2 dimensions (after Rodrigucs & Thomas 1995). Lower grid 
represents i-th level. Upper grids represent i + 1-th level with 
one of the overlapping grids displaced by half a blocklcngth. The 
hatched region in the lower grid is averaged to make a block at 
the i + 1-th level and indicated as the region marked as a thick 
square in the upper grid. Blocks in the offset grid are averaged 
from the lower level in the same way. 
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Figure 2. The overlapping criterion. The thick squares show a 
halo and the dashed squares show an investigating region. (a)The 
block can collapse and is merged into the halo because the over- 
lapping region is larger than half of the lesser of the regions of 
the halo and the investigating region. (b)The block cannot col- 
lapse because the overlapping region is not larger than half of 
the lesser of the regions of the halo and the investigating region. 
(c)Thc linking condition. In the case shown here, the block can- 
not collapse because we would like to prevent the growth of the 
long filamentary structure (see text) .Note that it is equivalent to 
consider either regions or masses for the overlapping criterion. 



how the mass function is changed. The condition (c), the 
Unking condition, prevents the growth of filamentary struc- 
tures because galaxies that we observe are not filamentary. 
However, we should note that the filament structures of dark 
matter certainly appear as shown in numerical simulations. 
Then the criterion for linking is just an assumption. 

It must be noted here that in the MCM, because over- 
lapping grids are used, the mass spectrum is close to contin- 
uous, rather than restricted to powers of two as in the Block 
model. 



2.2 Block model 

In this subsection, the Block model is briefly reviewed ac- 
cording to Cole & Kaiser (1988). Now, we consider a large 
block with mass Mo(~ W^^Mq), divide this into two, and 
realise blocks with various mass scales by successively divid- 
ing smaller and smaller blocks. We obtain a set of blocks of 
discrete masses. Mi = Mo/2*(i = 1, 2, • • • , N) where is an 
integer and is 27 in our calculation as shown later. 

Here a{Mi) = Oi means the standard deviation of the 
density fluctuations smoothed on a mass scale Mi. When the 
power spectrum is written in the scale-free form P{k) oc fc", 
the standard deviation is given by 



47rl/ 
(27i-)3 

M 

M, 



f 

Jo 



Wi{k)P{k)k'^dk 



(3) 



where M« is the mass on which the variance is unity, and 
Wi{k) is the window function. We define the quantities E, 
as follows: 



(To, 



a. 



2 



> 1. 



(4) 
(5) 



We generate the random density fluctuations on each 
block as follows: First of all, a density contrast generated by 
the random Gaussian distribution with the variance Eo is 
assigned to the largest block. This corresponds to the den- 
sity contrast smoothed on the mass scale Afo. Then the pos- 
itive random variable generated by the random Gaussian 
distribution with the variance Ei is added to one of the two 
divided blocks with the mass Mi = Mo/2 and a negative 
one with the same absolute value of the amplitude is added 
to the other divided block. Repeating this procedure from 
i — until i — N,wc obtain a tree of density contrasts of the 
blocks with various mass scales. The condition of collapse for 
a block at the redshift z is that the density contrast of the 
block equals 5c = 1.69(1 -|- z) {z; redshift). In this way, we 
can construct the merger tree of the dark haloes. 



3 ANALYTIC APPROACH TO MASS 

FUNCTIONS 

In this section, we present the PS formula, the Jedamzik 
formula and also the mass functions derived by our formula 
(YNG formula) which explicitly includes the effects of the 
spatial correlations. 



3.1 Press-Schechter formula 

The probability of finding the region whose linear density 
contrast smoothed on the mass scale M, 5m, is greater than 
or equal to 5c is assumed to be expressed by the random 
Gaussian distribution given by 



f{>Sc,M) 



27va{M) J,^ 



exp 



2cr2(M) 



d5. 



(6) 



This probability corresponds to the ratio of the volume of 
the region above 5c in the density contrast smoothed on 
the mass scale M to the total volume (in a fair sample of 
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the Universe). Therefore, the difference between /(> 5c, M) 
S'lid /(> 5c, M + dM) represents the volume of the region 
for which 5m = 5c precisely. The density contrast of an 
isolated collapsed object is precisely equal to 5c because an 
object with 5 > 5c would be eventually counted as an object 
of larger mass scale. The volume of each object with mass 
scale M is M/p. Then we obtain the following relation, 



Mn{M) 



dM 



df{>5c,M) 
dM 



dM, 



(7) 



where n{M) means the number density of objects with mass 
M, that is, the mass function. However, the underdense re- 
gions are not considered in the above equation. Hence, Press 
and Schechter simply multiply the number density by a fac- 
tor of 2, 



Ml<Ml,M^-2'Ji^^dM 
p dM 



(8) 



This factor of 2 has long been noted as a weak point in the 
PS formula (the so-called 'cloud-in-cloud' problem). Peacock 
& Heavens (1990) and Bond et al. (1991) proposed a solution 
to this problem by taking account of the probability P^p that 
subsequent filtering an larger scales might result in having 
(5 > (5c at some point, even when at smaller filters, 5 < 5c 
at the same point. They found that the factor of 2 in the 
PS formula could be correct only by using the sharp fc-space 
filter given by 



k < kc{Mi] 
k > kc{Mi] 



(9) 



3.2 Jedamzik formula 

Jedamzik (1995) proposed another approach to the cloud- 
in-cloud problem. 

Now, we consider the regions whose smoothed linear 
density contrasts on the mass scale M\ are above 5c- Each 
region must be included in an isolated collapsed object with 
mass M2 > Ml. Therefore, we obtain the following equation, 



/(><5c,Mi) 



P(Mi, M2)—n(M2)dM2, 

Ml P 



(10) 



where P{Mi, M2) means the conditional probability of find- 
ing a region of mass scale Mi in which 5mi is greater than 
or equal to 5c, provided it is included in an isolated over- 
dense region of mass scale M2. By 'the Jedamzik formula' 
we mean the procedure in which the mass functions are es- 
timated by solving eq.(lO). If P{Mi, M2) is given by the 
conditional probability p{5]\f^ > 5c\5m2 = Sc), P{M\, M2) is 
written as follows by using Bayes' Theorem, 

P(Mi,M2) = p(<5a/i > 5c\5m., = Sc) 
_ p{5mi > 5c,5m2 ~ Sc) 



p{5m2 = Sc) 

1 r r 1 {Sm, - Sc)^ 



27rcrsub Js^ 



2 '^Lb 



1 

2' 

'Tsub = a^(Mi)-a^(M2) 



dSMi 

(11) 
(12) 



can obtain the PS formula, naturally including the factor of 
2 as can be seen from eqs.(lO) and (pT|). 

However, it is insufficient for more realistic estimation of 
the mass function to use eq.(pT|) because it is necessary to 
consider the spatial correlation of the density fluctuations 
due to the finite size of the objects. Therefore, we must 
consider the probability P{r, M\, M2) of finding 5m-^ > 5c 
at a distance r from the centre of an isolated object of mass 
scale M2. Then we can get the probability P[M\, M2) by 
spatially averaging P{r, M\, M2)- 

Because we believe that the isolated collapsed objects 
are formed around density peaks, the constraints to obtain 
the above probability P{r, Mi, M2) are given as follows: 

(i) The linear density contrast, 5m2, of the larger mass 
scale M2, should be equal to 5c = 1.69 at the centre of the 
object. 

(ii) Objects of the mass scale M2 must contain a maxi- 
mum peak of the density field, i.e., the first derivative of the 
density contrast V5m2 must be equal to and each diago- 
nal component of the diagonalized Hessian matrix ^ of the 
second derivatives must be less than at the centre of the 
object. 

(iii) The density contrast of the smaller mass scale Mi(< 
M2) which collapsed and is included in an object of mass 
scale M2 must satisfy the condition 5mi > Sc at distance r 
from the centre of the larger object. 

The probability which we get from the above conditions is: 
P(r,Mi,M2|peak) = P{SM^ > ^d^Afa , peak) 



27r(l - e2 _ ^2 _ ^2 2e^i7) 

/o°° fi^) J^^ dui exp(- 

/;°dx/(a;)exp(-%) 



(13) 



Since the detailed derivation of the above probability and 
explanations of the notation are complicated, we omit them 
here (see YNG) . By spatially averaging eq. (M± in the region 
of M2, we obtain P(Mi,M2) as follows: 



P(Mi,Af2 



/g^' P(r, Ml , Af2 |peak)47rr2dr 
4nr^dr 



(14) 



where we use the sharp fc-space filter (see YNG). Thus we 



where R2 means the radius of the region M2, R2 ~ 
{3M2/4:np)^^^^\ The cumulative multiplicity functions, 
which will be defined in the next section, are estimated by 
using eq.(^3|) and are shown by the short-dashed lines in 
Fig. 3 in the cases of the power spectrum P{k) oc fc" with 
n — and —2, respectively. We call this formula the YNG 
formula hereafter. 



4 RESULTS 

4.1 Cumulative multiplicity functions 

We calculate the cumulative multiplicity function P{> M), 
which is defined as the mass fraction of objects whose 
mass is larger than a mass M to the total mass (in a fair 
sample of the Universe), by following the various methods 
which are mentioned in §§2 and 3. The multiplicity function 
P{M)d\nM (which we define in a logarithmic interval of 
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n=0 




.2 .4 .6 

log(M/M.) 

n=-2 



A 



.2 .4 .6 .8 1 

log(M/M.) 

Figure 3. Cumulative multiplicity function: (a)spectral index 
n=0, (b)n=-2. The solid lines, the short-dashed lines, the dotted 
lines and the long-dashed lines show the cumulative multiplicity 
functions given by the MCM averaged over four realisations, the 
YNG formula, the PS formula and the Block model averaged over 
five realisations, respectively. 



mass) and the cumulative multiplicity function are related 
to the mass function n{M)dM as follows: 



P{M)dlnM 



dp 
P 



IVPn(M) 



dlnM, 



P(>M) = r P{M)^. 

J M 



(15) 
(16) 



We assume that the linear density fluctuations obey a ran- 
dom Gaussian distribution with a power-law power spec- 
trum P{k) oc fe", where n = and —2. We consider only 
the Einstein-de Sitter universe (f2 = 1,A = 0) in this pa- 
per. In the MCM, we take the box size L = 128. M, 
is deflned as cr(Mt) = 1 (see eq.(3)) and here the mass 
in the block with the mass scale M, is assigned to eight 
cells. In the Block model, the largest box size is assigned 
to the mass scales Mo = W M-^ in the case of n = and 
lO^M* in the case of n = —2. We consider block sizes with 
Mi = Mo/2'(i = 1,2, = 27). 

In Figs. 3(a) and (b), we show the cumulative multiplic- 
ity functions derived by the MCM, the YNG formula, the 
PS formula and the Block model for the case of n = and 
n = 2, respectively. 

The multiplicity functions of the MCM are shown only 
in the mass range (0 < log(M/M,) ^ 1) which corresponds 
to the range from eight cells to about 100 cells. The reason is 
as follows: On the smallest mass scale (one cell scale) in the 
numerical calculations, the power of the density fluctuations 



on the one cell scale cannot be correctly produced in agree- 
ment with the theoretically estimated power of the scale-free 
mass spectrum a{M) oc M"^^"*""^/® due to the finiteness of 
the cell size in the numerical calculations. On the larger 
mass scales, larger haloes than about 100 cells have various 
shapes, rather than the spherical shape, so the assumption 
of spherically symmetric collapse fails. Hence, the cumula- 
tive multiplicity function given by different method on these 
scales cannot be directly computed without clarifying the 
identification of the isolated haloes. Moreover, since there 
are few haloes on the scales larger than about 100 cells in 
the numerical calculations, the error in the number of haloes 
increases. Therefore, wc show here the cumulative multiplic- 
ity function only on the mass range from eight cells to about 
100 cells (0 < log(M/M,) <, 1). 

It is found that the multiplicity functions given by the 
MCM are well flt by those given by the YNG formula. In 
the case of n = 0, the agreement is good while in the case 
of n = —2, there is a little deviation on the mass scale 
log(M/M,) ^ 0.8 because of the numerical errors mentioned 
above. It is also found that the PS formula and the Block 
model are also in agreement with each other in the case of 
n = —2 rather than in the case of n — 0. Note that on the 
galaxy scales the spectral index n is nearly equal to —2 in 
the CDM model. Then in the CDM model the both func- 
tions are consistent with each other (Cole & Kaiser 1988). 
It must be noticed that the agreement of the Block model 
on the PS formula depends on the power spectrum. On the 
other hand, the function of the MCM deviates from those 
of the Block model and the PS formula in which the spar 
tial correlations of the density fluctuations are not explicitly 
taken into account. From these results, we can conclude that 
the MCM naturally and correctly takes into consideration 
the spatial correlations and the deviation of the multiplicity 
function of the MCM from those of the Block model and the 
PS formula results from the effect of the spatial correlations 
of the density fluctuations. Furthermore, we will show in the 
next subsection that the multiplicity function given by the 
Block model can be reproduced by using the Jedamzik for- 
mula without consideration of the spatial correlations and 
thereby demonstrate that the difference between the mul- 
tiplicity functions given by the Block model and the MCM 
results from the effect of the spatial correlations. 



4.2 Block model 

Here we analytically reproduce the multiplicity function 
given by the Block model by using the Jedamzik formula. 

Since we consider the discrete mass of the blocks {Mi = 
Mo/2') in the Block model, the density contrast in the blocks 
which are identifled as isolated collapsed haloes is gener- 
ally greater than S^. We cannot recognise the just collapsed 
halo whose density contrast is precisely Sc if we follow the 
procedure of the Block model. Therefore, the conditional 
probability P(Mi,M2) in this case must be approximately 
expressed as: 

P(Ml,M2) = p(Sm, > S,:\Sm2 > S,:) 

p{Sm2 > <5c) 
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Figure 4. Block Model: (a)n=0, (b)n=-2. Crosses, solid squares 
and open squares show the multiplicity functions averaged over 
five realisations of the Block model with the largest box size A/q = 
10^M,,10^M, and lOM* in the case of n = and lO'^Af* , 10^ Af, 
and lO^Af* in the case of n = —2, respectively. The solid lines 
and the dotted lines show the predictions by using the Jedamzik 
formula and the PS formula, respectively. 



^ SZdu.^j----v[-%^-i} ^^^^ 

v/2.(l-e'2)/-d^2exp{-|} 

It must be noticed here that we must consider the case of 
^ instead of the case of (5 a/2 = which appeared 
in eq.(ll). In this case, of course, the spatial correlation is 
not taken into consideration. The detailed derivation of the 
above probability and notation are explained in Appendix 
A. Inserting the above conditional probability, eq.(17), into 
eq.(lO), we can estimate the multiplicity function in this 
case. 

In Fig. 4, we show the multiplicity functions given by 
the above procedure (solid line) and the Block model (cross, 
solid square and open square). The cross, solid square and 
open marks in the Block model mean that the masses of the 
largest block Mo equal 10^M*,10^Af, and lOA'/*, respec- 
tively. Here in the Block model, the multiplicity function 
over a large dynamic range in mass can be produced by 
combining the functions derived by the Block model with 
Mo = l<fM,,\Q^M-^ and lOAf. in the case of n = and 
IffM,, 10^ Af, and 10^ A/, in the case of n = -2 and block 
sizes with Mi = A/'o/2'(j = 1, 2, • ■ • 27) in each case. We find 
that the multiplicity functions given by the Block model and 
the analytical formula derived from the Jedamzik formula 
fit each other well. We believe that the difference between 



the cumulative functions given by the Block model and the 
MCM mainly results from the effect of the spatial correla- 
tions while the difference between the functions given by the 
Block model and the PS formula results from the difference 
in the identification of isolated collapsed objects. 

4.3 Overlapping effect 

In taking into consideration the finite size of the haloes, it 
appears necessary to consider the serious problem of the 
spatial overlapping of the dark haloes. For the overlapping 
criterion we must consider how we can identify the number 
and sizes of the haloes when some haloes overlap. In this 
subsection, we show the effect of the overlapping criterion 
(§2.2,(b)) on the multiplicity function. The overlapping cri- 
terion adopted in RT is that the investigating region can 
collapse when the overlapping region has a larger mass than 
half of the lesser of the masses of the halo and the investi- 
gating region. Here we quantify the overlapping criterion by 
defining a parameter, x, as the ratio of the mass of the over- 
lapping region to the lesser of the masses of the halo and 
the investigating region. Then, RT's criterion corresponds 
to a; = 1/2. By changing the value of x we show the effect 
of overlapping on the multiplicity function. 

In Fig. 5, we show the multiplicity functions on scales 
above eight cells in the cases that x — 1/2, 1/4 and 1/8. As 
the value of x decreases, the multiplicity function increases 
on the larger mass scales because it becomes easy for larger 
blocks to collapse. 

We find that the multiplicity function on the larger 
mass scales changes according to the value of x. Note that 
this change is shown quantitatively in Fig. 5 only on scales 
smaller than about 100 cells because of the uncertainties 
due to the numerical errors on scales larger than 100 cells 
which are mentioned in §4.1. We consider the trend of the 
multiplicity function on large mass scales to increase with 
decreasing x on scales larger than 100 cells to also be signif- 
icant. 

So, the overlapping effect is a serious problem for con- 
structing merger trees of dark haloes, particularly on mass 
scales larger than about 100 cells. 



5 CONCLUSIONS & DISCUSSIONS 

In this paper, we have shown that the multiplicity functions 
given by the MCM are consistent with those given by the 
YNG formula. However, the functions of the MCM do not 
fit those given by the PS formula and the Block model. This 
fact means that the MCM includes the information of the 
spatial correlations in the density fluctuation field, but the 
PS formula and the Block model do not include it. Thus we 
believe that the effects of the spatial correlations affect the 
merger trees as well as the mass functions and it is important 
to take into account the spatial correlations in the merger 
tree models. 

Here it must be noted that the reason why the multiplic- 
ity functions given by the MCM are well fit by those given by 
the YNG formula in the case that the overlapping criterion 
has the value of a; = 1/2 is as follows. When obtaining the 
multiplicity function from eqs.(^3|) and ([l^), we integrate 
the P{r,Mi, M2) with respect to r from to R2 in eq.(|l^. 



© 0000 RAS, MNRAS 000, 000-000 



8 M. Nagashima & N. Gouda 



n=0 




Figure 5. Effects of overlapping: (a)n=0, (b)n=-2. Tfie solid 
lines, dotted lines and dashed lines show the multiplicity func- 
tions averaged over four realisations given by the MCM with the 
overlapping criterion defined hy x = 1/2,1/4 and 1/8, respec- 
tively. 



This corresponds to the case that x — 1/2. Of course, the 
adoption of the integration interval [0, R2] is, strictly speak- 
ing, only justified for spherically symmetric collapse. In fact, 
however, the effects of non-spherical collapse or tidal inter- 
action among haloes may be important. We hypothesise that 
changing the value of x might be enough to account for these 
effects. 

We also analytically produce multiplicity functions 
which fit those given by the Block model, by using the 
Jedamzik formula. When we reproduce the multiplicity func- 
tio ns g iven by the Block model, the conditional probability 
eq.(|r^ is adopted instead of eq.(^l|) in the PS formula. The 
difference from the PS formalism is that the density con- 
trast of the isolated object is not only equal to Sc but also 
greater than Sc. This difference results from the discreteness 
of block mass in the Block model. It must be noticed that 
the conditional probability eq.(^) does not include the ef- 
fect of the spatial correlations of the density fluctuations. 
We find that the multiplicity functions reproduced in this 
way are consistent with those given by the Block model. 
This result shows that the Block model does not include the 
spatial correlations of the density fluctuations and so the 
multiplicity function given by the Block model is different 
from those given by the MCM. 

Moreover, we have shown that when the overlapping 
criterion is changed, the multiplicity function is affected, 
especially on larger mass scales. As the value of the overlap- 
ping parameter, x (defined as the ratio of the mass of the 



overlapping region to the lesser of the masses of the halo 
and the investigating region) decreases, the number density 
of high mass haloes increases because it becomes easy for 
larger blocks to collapse. 

We believe that the MCM is a powerful tool for con- 
structing merger trees of dark haloes because it includes 
the information of spatial correlations naturally and cor- 
rectly as shown in this paper. In the MCM, the number of 
high mass haloes increases, relative to the results of previous 
works, e.g., Kauffman et al. (1993) and Cole et al. (1994), 
so the number of giant, red elliptical galaxies may increase. 
In the MCM, however, a serious problem appears, i.e., it 
is found that in taking into consideration the overlapping 
of dark haloes of non-zero size, the multiplicity function is 
sensitive to the choice of the overlapping criterion. There- 
fore, we must consider a more realistic overlapping criterion 
for constructing merger trees. We do not know a priori the 
value of the overlapping parameter. Furthermore the value 
of the overlapping parameter may not be constant as time 
passes. It might be a function of halo mass, of separation 
between a halo and a block, or some other parameter. It 
is a very difficult problem to identify the mass and size of 
the collapsed halo when several haloes are overlapped. How- 
ever, it is very important to precisely identify the collapsed 
haloes for constructing the merger trees of haloes with good 
accuracy. As mentioned in §1, the merger tree of the dark 
haloes greatly affects the various important processes for the 
formation and evolution of galaxies, such as the thermal his- 
tory of gases, merging between galaxies and so on. Hence, it 
is necessary and important to understand how the overlap- 
ping criterion influences these processes. Furthermore, we 
assume spherically collapse of the dark haloes. But, in fact, 
there is the possibility that the collapsed objects have flla- 
ment or sheet structures. We must also deal with these cases 
in order to correctly estimate the merger trees. We are now 
investigating the merger tree models in which the effect of 
non-spherical collapse is included and the identification of 
the collapsed haloes is well defined. 
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APPENDIX A: MASS FUNCTION OF BLOCK 
MODEL 

As mentioned in §3, we show the analytic formula which 
reproduces the mass functions given by the Block model. 

The analytic formula is approximately expressed as fol- 
lows. In YNG and PS, it is assumed that the density contrast 
5m2 of the isolated object which has collapsed with mass A/2 
is just 5c = 1.69 (see eq.(13)). In the Block model, however, 
this condition should be changed to SAt^ > Sc because in 
general we can not see the block which has just collapsed 
and its density contrast is just Sc due to the discreteness 
of the mass scales of the blocks with A/i+i = Mi/2 in the 
Block model. 

So the probability should be changed as follows; 



P(Mi,M2) 



v(Smi > 5c\5m2 > 5c) 
p{Smi > Sc, Sm2 > Sc) 



(Al) 



p{Sm2 > Sc 

In orde r to estimate the probability on the right hand side of 
eq.(Al), we have to consider the two- variables Gaussian dis- 
tribution function. The probability of m-variables Gaussian 
is generally (see BBKS) 

-Q/2) 



p(V)dV = 



^(27r)'"det(M) 



dV, 



where 

Q - 

V 

Mii ■- 



VM^^V"^, 

{yi,y2, ■ 



iivi - {yi))iy3 - (w>)> 



(A2) 



(A3) 
(A4) 
(A5) 



and yi{i — 1,2, ... ,m) is the Gaussian random variables. 
We use the angle brackets {) as the ensemble average of the 
universe, but in practice, assuming homogeneity and ergod- 
icity in space, we can take it as the spatial average. So (yi) 



corresponds to the spatial average of yi, and Mij is covari- 
ance between yi and j/-, . In this case, we may consider the 
two variables, Sm^ and Sm2, so the covariance matrix M is 
written as follows: 



jyj ^ [ (S'lh ) (Smi ) 

{SM1SM2) (Sli^) 



(A6) 



Note that {i5mi) s^id (^Mi) ^-r^ variances, a^{Mi) and 
a^{M2), respectively, and {Shi-^Sm2) is a cross correlation 
at the same point. 

Here, we normalize the density contrast and the cross 
correlation, 

al = {SmiSm2). (A7) 



a {Mi] 



a{Mi)a{M2)' 



Using the above notation, we obtain the two variables Gaus- 
sian distribution function, 



1 



2-kVI^ 
So we obtain 



: exp 



[vi - eV2)'^ _ V2_ 
2(1- £2) 2 



dv\ dv2 ■ 



(A8) 



P('5aIi > Sc,Sm2 > Sc) =P{^1 > Vlc,V2 > !^2c) 
1 



27r^ 



- 

00 /• oc 



T X / / cxp 



[vi - eu2 



2(1- £2) 2 



^ ) duidu2{A9) 



where fie is Sc/(j{Mi). Furthermore, dividing the above 
equation by the integration of one variable Gaussian dis- 
tribution. 



P{Sm2 > Sc) = P{f2 > f2c) 
1 



2. ./.^/'Pi"?''^'''' 



(AlO) 



we obtain the P(Mi, M2). 

Next, we estimate the normalized cross correlation func- 
tion e. In eqs.(^ and (^, the variance ai with the mass Mi 
is given by 



AtvV 



Wi{k)P{k)k^dk, 



(All) 



and the variance of the additional density contrast E is given 
by 



Si — o\ — af_i, i > 1. 



(A12) 



Because we give the additional density contrast indepen- 
dently of another density contrast with the other mass scale, 
the variance of the additional density contrast must depend 
on the power spectrum only at the interval of the wavenum- 
ber [fci+i,fei], corresponding to the mass scales [Mi, Mi+\\. 
So we should adopt the sharp- space filter as the window 
function Wi{k) defined by, 



W,{k) = 



1, 



k < kc{Mi) 



(A13) 



0, k> kc{Mi), 
In this case, the variance and the cross correlation are 



(AM) 



© 0000 RAS, MNRAS 000, 000-000 



10 M. Nagashima & N. Gouda 



^ft = / Wi{k)W2{k)P{k)k''dk 



(27r)3 
471-y 



fc2 

P{k)k'^dk 



(27r)3 

= ol (A15) 

where i = 1,2 and A/i < A/2, i-c, k\ > k2. So the normal- 
ized cross correlation function e is 

e=^^^. (A16) 

0"1(T2 (Tl 



Substituting the above eq.(A16) into eq.(A9), we obtain 



oc /• oc 



Ac "'ic \ sub / 

where 

'^sub = '^i~'^2- (A18) 
Therefore the conditional probability P(A/i, ^^2) is 

P{Mi,M2)^N-\ i 

27ro-2crgyjj 

Sc JSc \ ^''^sub 



c^P - 02 - TT^ '^^lO^'^a, (A19) 



1 /-^ / SI 



where 



Note that this probability is the same as eq.(8) in Jedamzik 
(1995). 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M^rjK style file. 
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